Penetration depth, multiband superconductivity, and absence of 
muon-induced perturbation in superconducting PrOs 4 Sbi 2 



Lei Shu*, D. E. MacLaughlin, and W. P. Beyermann 
Department of Physics, University of California, Riverside, California 92521 

R. H. Heffner 

Los Alamos National Laboratory, Los Alamos, New Mexico 87545 

G. D. Morris 

TRIUMF, 4 004 Wesbrook Mall, Vancouver, B.C., Canada V6T 2A3 

0.0. BernaF 
Department of Physics and Astronomy, 
California State University, Los Angeles, California 90032 

F. D. Callaghan and J. E. Sonier 
Department of Physics, Simon Fraser University, Burnaby, B.C., Canada V5A 1S6 

W. M. Yuhasz, N. A. Frederick, and M. B. Maple 

Department of Physics and Institute for Pure and Applied Physical Sciences, 
University of California, San Diego, La Jolla, California 92093 
(Dated: April 21, 2009) 



* Current address: Department of Physics and Institute for Pure and Applied Physical Sciences, University 

of California, San Diego, La Jolla, California 92093 
' Current address: Division of Materials Research, National Science Foundation, Arlington, VA 22230 



1 



Abstract 

Transverse-field muon spin rotation (TF-//SR) experiments in the heavy-fermion superconductor 
PrOs4Sbi2 (T c = 1.85 K) suggest that the superconducting penetration depth A(T) is temperature- 
independent at low temperatures, consistent with a gapped quasiparticle excitation spectrum. In 
contrast, radiofrequency (rf) inductive measurements yield a stronger temperature dependence 
of A(T), indicative of point nodes in the gap. Muon Knight shift measurements in the normal 
state of PrOs4Sbi2 suggest that the perturbing effect of the muon charge on the neighboring Pr 3+ 
crystalline electric field is negligibly small, and therefore is unlikely to cause the difference between 
the TF-;uSR and rf results. The discrepancy appears to be related to multiband superconductivity 
in PrOs 4 Sbi 2 . 

PACS numbers: 71.27.+a, 71.70.Ch, 74.70.Tx, 75.10.Dg, 76.75.+i 
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I. INTRODUCTION 



In the phenomenon of multiband superconductivity (MBSC), which was first treated 
theoretically in 1959 using BCS theory^ distinct energy gaps open up on different sheets of 
the Fermi surface in the superconducting state. One of the well-studied multiband super- 
conductors is the binary intermetallic compound MgB2, 2,3 ' 4 which has two bands and two 
superconducting gaps. Based on a quasiclassical theory, microscopic calculations of the elec- 
tronic structure in the vortex state have been made in a two-band superconducting model.- 1 ^ 
These suggest that at low applied fields the dominant contribution to the change of the total 
density of states (DOS) comes from the small-gap band, which gives rise to spatially ex- 
tended quasiparticles (QPs). At high fields these loosely bound states become delocalized, 
with the vortex core size determined by the more localized states associated with the larger 
gap. The presence of delocalized QPs modifies the DOS, transport properties associated 
with the quasiparticles, and the spatial field distribution around a vortex.- 

Recently strong evidence for MBSC has been found in the filled skutterudite PrOs4Sbi2 
from thermal transport measurements in the vortex state.- 1 ^ PrOs4Sb 12 has attracted much 
attention since its superconductivity was discovered in 2002.-- It is the only known Pr- 
based heavy-fermion superconductor (T c = 1.85 K) and exhibits a number of extraordinary 
properties.— The Pr 3+ ground state is a nonmagnetic Ti singlet, which is separated from a 
I4 first excited state (tetrahedral notation^) by only ~10 K.— iMii£ii& A novel ordered phase 
appears at high fields and low temperatures. There is evidence that time-reversal symmetry 
is broken in the superconducting state.— Radiofrequency (rf ) inductive measurements of the 
magnetic penetration depth A in the Meissner state,— thermal conductivity measurements 
in a rotated magnetic field,™ and flux-line lattice distortion 2 ^ all suggest nodes in the su- 
perconducting gap. A double superconducting transition has been observed in specific heat 
measurements.— i2ii22i22 However, recent specific heat and heat transport measurements^ on 
a highly homogeneous single crystal show only one transition peak in the specific heat and 
a fully-gapped Fermi surface. The latter result is corroborated by angle- dependent specific 
heat measurements on single crystals.— Previous transverse-field muon spin rotation (TF- 
/iSR) measurements of A in the vortex state of PrOs4Sbi2 found evidence for a BCS-like 
activated dependence at low temperature—? suggesting the absence of gap nodes. Thus 
there are a number of open questions concerning the superconducting order parameter in 
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PrOs4Sbi2- 

The present article reports a detailed TF-/1SR study of PrOs4Sbi2- TF-/zSR experiments^ 
have proved invaluable in characterizing both the superconducting and normal states of 
superconducting materials;— 1 ^ in particular, a TF-/zSR studyi 1 of superconducting NbSe2 
clearly revealed effects of MBSC on the vortex-state field distribution in this compound. 
We compare TF-/zSR and rf inductive penetration depth measurements in superconducting 
PrOs4Sbi2, and discuss the previously- reported^ 2 - discrepancy between these measurements 
in terms of extreme MBSC in PrOs4Sbi2- Preliminary reports of parts of this work have 
been published.— 1 ^ 

We have also carried out TF-/iSR measurements of the Knight shift in the normal state 
of PrOs 4 Sb 12 , which suggest that the perturbing effect of the positive-muon charge on 
the neighboring Pr 3+ crystalline electric fields (CEF) is negligibly small. This indicates that 
the muon charge is unlikely to be the source of the discrepancy between the TF-/iSR and 
rf results. A model calculation for the perturbed CEF energy levels, based on the approach 
of Kaplan, Schenck, and co-workers,— is&2£i2Z is described. Higemoto et al. 38 have reported 
a TF-/zSR study of PrOs4Sbi2, with emphasis on the Knight shift in the superconducting 
state. 

The remainder of this introduction contains three brief pedagogical sections: a descrip- 
tion of the elements of the TF-/zSR technique used in this study (Sec. II Aj) . a review of 
the important features of the /i + Knight shift (Sec. II Bp . and an introduction to the CEF 
model calculation (Sec. II CI) . After a description of the experimental procedure (Sec. [TTJ) , in 
Sec. II I II we describe our experimental results in PrOs4Sbi2, which include the temperature 
dependence of the TF-/1SR relaxation rate and penetration depth, magnetic susceptibility 
data, and normal-state /i + Knight shift measurements. The implications of these results for 
the nature of both the superconducting and the normal state are discussed in Sec. [TV] We 
summarize our results in Sec. [V] 

A. Transverse-field muon spin rotation in the vortex state 

In the TF-/xSR technique 25 spin-polarized positive muons (/i + ) are implanted in the sam- 
ple and precess in a local magnetic field applied perpendicular to the initial fi + polarization. 
This precession is detected using the asymmetry of the fi + beta decay (the decay positron 
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is emitted preferentially in the direction of the /i + spin). The distribution of fi + precession 
frequencies directly reflects the distribution of magnetic fields in the sample. Thus TF-/iSR 
can be used to measure the magnetic field distribution of the vortex lattice in a type-II 
superconductor and the local magnetic susceptibility in the normal state. 

In the vortex lattice each vortex possesses a normal-state-like core with radius of order 
of the superconducting coherence length £, surrounded by a shielding supercurrent. These 
supercurrents give rise to an inhomogeneous magnetic field that is periodic in the vortex 



lattice. An analytical Ginzburg-Landau (GL) model for the spatial field profile B(r) (Ref . l39l) 
yields 

B(r) = ft<l-*)£f^<H>*, (1) 

K 

where the K are the reciprocal lattice vectors of the vortex unit cell, K\{u) is a modified 
Bessel function, b = B / ' B c2 is the reduced field, and 

u 2 = 2£ 2 K 2 (l+b 4 )[l-2b(l-b) 2 }, (2) 

with £ the GL coherence length. This result is valid in the limit X 2 K^ in ^> 1, where K min is 
the smallest nonzero reciprocal lattice vector. 

TF-/iSR is a sensitive probe of this field distribution.— The functional form of the /i + spin 
relaxation, which depends on the field distribution, is fit to the functional form G(t) cos(co> M t+ 
<p) , where the frequency uj^ and phase <f> describe the average fi + precession and the relaxation 
function G(t) describes the loss of phase coherence due to the distribution of precession 
frequencies. The relaxation rate associated with G(t) is a measure of the rms width SB of 
this distribution. The expression- 

5B 2 (T) = 0.00371$^A" 4 (r), (3) 

where <J>o is the flux quantum, then gives a rough estimate of A for a triangular vortex lattice 
in the London limit (A>^). More accurate relations are available,— £&£k but in the present 
case lead to negligible corrections to Eq. (j3J). 

To estimate the field distribution width without fitting a theoretical model, a Gaussian 
distribution of local fields is often assumed. Then the time dependence of the fi + spin 
polarization is proportional to exp(— ^a 2 t 2 ), where a = j^SB is the /i + relaxation rate; 7 M is 
the fi + gyromagnetic ratio. It has been pointed out^ that this procedure is approximate at 
best, because the field distribution is not expected to be Gaussian. In a more microscopic 



approach, a lineshape analysis program has been written to fit the GL model to TF-/iSR 
data.— ^ The GL model is preferred over the London model because it self-consistently 
accounts for the vortex cores. In this fitting program the fit function is the Fourier transform 
of the magnetic-field distribution function n(B') = J d 2 r 5[B f — B(r)], which can be obtained 
from Eq. ([T]). These fits utilize the entire form of the field distribution, not just its second 
moment. They yield an effective '/zSR' penetration depth A m sr, which becomes the true 
penetration depth only after extrapolation to B = 0.— 



Results from such fits in PrOs 4 Sbi 2 (Ref. l32[) are described below in Sec. IIII Al 



B. Muon Knight shift 

In TF-/xSR the total field at the /i + site is given by the sum of the applied field H, the 
internal field induced by the applied field, and the demagnetization and Lorentz fields™, 
The relative /i + frequency shift 

k; = ^-i, (4) 

where u; re f = 7^ H is the /i + frequency in "free space" (i.e., no condensed matter effects), 
must be corrected for the contribution of the demagnetization and Lorentz fields to 
obtain the /i + Knight shift^ 3 - 

k, = k;- K Dh , (5) 

which contains the relevant information about the local magnetic susceptibility. 

In a paramagnetic metal, originates from hyperfine fields produced by the field- 
induced polarization of conduction electrons and localized electronic moments. The con- 
tribution from the conduction electrons is temperature independent and is usually very 
small, of the order of 100 ppm. 43 The local moments, in the present work the Pr 3+ ions, 
contribute to via two coupling mechanisms: (1) the dipolar interaction between the local 
moments and the which may be described as a dipolar field at the /i + interstitial site, 
and (2) the indirect RKKY interaction,— in which an additional spin polarization of the 
conduction electrons due to the local moments produces a hyperfine contact field at the 
interstitial fi + . Both contributions are proportional to the local-moment susceptibility: 

K; = (A^ + A con ) Xl (i = x,y,z) (6) 
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in the principal- axis coordinate system of the dipolar tensor Adi P , where the A^ ip are the 
diagonal elements of Adi P and A con is the contact hyperfine coupling constant, assumed 
isotropic. For a cubic lattice (x x = Xy — Xz), only the hyperfine contact field at the 
interstitial /i + site contributes to the average shift = | ^\ since the contribution 
from the dipole-dipole interaction vanishes A 1 ^ = 0) . 

The local susceptibility is sensitive to change in the Pr 3+ CEF, and the /i + charge may 
induce such a change. If so, the modified local susceptibility will be reflected in a breakdown 
of the proportionality of the yU + Knight shift to the measured bulk susceptibility ^ bulk since 
in this case x-i x bulk - This effect has been studied in detail by Kaplan, Schenck, and 
co-workers^ 1 ^ 1 ^ 1 ^ in Pr-based compounds with singlet Pr 3+ CEF ground states. As noted 
above and discussed in Sec. II C\ our results in PrOs4Sbi2 show little if any such effect. 

C. CEF model calculation 

The crystal structure of PrOs4Sbi2 belongs to the Im3 space group, with Pr 3+ ions at 
the points of a bcc unit cell. As an initial approximation we assume the CEF Hamiltonian 
Hcef for the Pr 3+ ions has Oh (cubic) point group symmetry. Then 

Hcef = BlO\ + hBlOl + B° 6 0° 6 - 215 6 °0 6 4 , (7) 

where 0°, Of, 0$, and 0\ are the Stevens operators for a given angular momentum quantum 
number J, and the Bs are parameters usually determined from experiment.— 7Ycef splits 
the Pr 3+ J = 4 Hund's-rule multiplet into a Ti singlet, a T3 doublet, and two triplets and 
r 5 ). Pr 3+ ions in PrOs4Sbi2 actually have (tetrahedral) point group symmetry, however,— 
for which an additional sixth-order term appears in the CEF Hamiltonian [Eq. (17!)]. This 
term mixes the T 4 and T 5 triplet wave functions with each other— but has a relatively small 
effect on the CEF energies.— If it is predominant, a T 4 or T 5 triplet will be the ground state, 
inconsistent with the singlet ground state in nonmagnetic PrOs4Sbi2. However, if this term 
is small, the change in the physical properties in zero magnetic field can be approximated 
by changing the parameters in the cubic CEF Hamiltonian [Eq. (J7j)] slightly. Thus we use 
the cubic CEF Hamiltonian for simplicity. 

In the presence of an external magnetic field, the Zeeman interaction mixes and splits 
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the CEF energy levels. The magnetic susceptibility is given by 37 



XCEF — 



y 

Z— in 


~(E n iy ) 2 /kT-2E n 2) 


exp ( 


-E n 0) /kT) 


£ n exp (-E^/kT) 





(8) 



where the E n are the unperturbed cubic CEF levels, E n = /J,Bg{4>n\J\(j>n) with the CEF 
wave functions and g the Lande g-factor, and 

K'-Z^^g—^ — -^y- (9) 

i , -C/n — / 

n =£n n 

In the molecular-field approximation, the measured magnetic susceptibility \ is given by 
X — Xcef/(1 — £xcef), where £ is the molecular field parameter that describes exchange 
interactions between Pr 3+ ions. 45 

Equation (jSJ) shows that the observed susceptibility is directly related to the local CEF 
energy levels. As noted above, if we take into account possible muon-induced change of the 
CEF, we must consider a modified CEF Hamiltonian to calculate the resultant change of the 
Pr 3+ local susceptibility— li&i&Sl In PrOs4Sbi2 the Pr 3+ first excited state is separated from 
the singlet ground state by only ~10 K.— liiiiLiS It is particularly important to determine 
possible muon-induced changes in the local CEF energy levels, since CEF excitations are 
central to a number of theories of superconductivity in PrOs4Sbi2-— ^ Changes in CEF 
splitting might therefore affect local superconducting properties such as the vortex-state 
super current. 

II. EXPERIMENTS 

TF-//SR experiments were carried out at the M15 and M20 channels at TRIUMF, Van- 
couver, Canada, on a mosaic of oriented PrOs4Sbi2 crystals. The crystals were mounted 
on a thin GaAs backing, which rapidly depolarizes muons in transverse field and minimizes 
any spurious signal from muons that do not stop in the sample. TF-/xSR asymmetry data3^ 
were taken for temperatures in the range 0.02-250 K and hqH between 10 mT and 1.0 T 
applied parallel to the (100) axes of the crystals. For the /i + Knight shift measurements, the 
applied field was determined by measuring the precession frequency of muons that stopped 
in a small piece of silver foil included with the sample. 
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For the superconducting-state TF-/xSR measurements the magnetic field H was applied 
in the normal state, and the sample was cooled in constant field to below T c . For PrOs4Sfc>i2 
H c i(0) < 45 Oe,— so that H was significantly larger than H cl and the sample was always in 
the vortex-lattice state. The demagnetizing field — 4ttDM, where D is the demagnetization 
coefficient and M < in the superconducting state, adds to the field B = H + 4n(l — D)M 
inside the sample, rendering it closer to the applied field. In Sect. IHI CI we estimate an 
effective D ~ 0.8, which yields an estimated variation of B with H c \ (i.e., with temperature) 
of only a few percent. The strong bulk vortex pinning in PrOs4Sbi2 (Ref. 1281) also suppresses 
changes in B. Thus modification of the vortex-state field distribution by flux expulsion, 
demagnetization effects, etc., seems unlikely to play an appreciable role in the experiments. 



III. EXPERIMENTAL RESULTS 



A. Superconducting-state TF-/iSR and magnetic penetration depth 

In this section we describe the temperature and field dependencies of the TF-/iSR re- 
laxation data in superconducting PrOs4Sbi2, and compare the superconducting penetration 
depth A(T) obtained from these data with inductive measurements^ in the Meissner state. 
We concentrate on A(T) at low temperatures, where power-law behavior is evidence for gap 
nodes and temperature-independent (or activated) behavior signals a fully-gapped Fermi 
surface. The behavior of the data at higher temperatures is more complicated and will 
not be considered in detail. Furthermore, as noted above in Sec. II A\ the effective TF-//SR 
penetration depth is expected to be field dependent, and approximates the true value only 
as H — > 0.— Thus we concentrate on the results at low temperatures and (for the TF-/xSR 
data) low fields. 

Figure [U compares A m sr(T) in PrOs4Sbi2 at low temperatures, obtained from GL model 
fits to TF-/zSR data taken at HqH = 10 mT, with the penetration depth increase AA r f(T) = 
A(T) — A(0) in the Meissner state obtained from rf inductive measurements.-^ The rf data 
do not determine A(0) separately; in Fig. [I] the quantity plotted (triangles) is A r f(T) = 
AA r f(T) + A a1 sr(0), which allows direct comparison of the temperature dependencies from 
the two experiments. 

It can be seen that X^sr(T) is nearly constant below ~0.6 K, indicative of a fully gapped 
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FIG. 1: (color online) Temperature dependence of penetration depth A in PrOs4St>i2 below ~1 K. 
Squares: A^sr from GL model fitting. Curve: fit to BCS low-temperature expression A(T) = 
A(0)[1 + y / 7rA/2Te- A / T ], A(0) = 0.3534(24) fj,m, 2A/k B T c = 4.9(1). Triangles: A rf from rf 
inductive measurements 1 ^ in the Meissner state (see text). Curve: fit to power law A + BT n , 
n = 2.05(3). Inset: Temperature dependence of A^sr up to T C M Squares: A^sr as in main panel. 
Circles: A^sr from Gaussian fits and Eq. (fTT|) . 



quasiparticle excitation spectrum. As previously noted^ the BCS low-temperature expres- 
sion A(T) = A(0)(1 + v/7rA/2Te- A / T ) (curve in Fig. [I]) gives a good fit to the TF-/iSR data 
for T < 0.5T C , suggesting that the energy gap is isotropic. The increase of A r f(T) (triangles) 
with increasing temperature is much stronger than the increase of A m sr(T) and follows a 



power law A + BT n , n « 2 (Ref. 1181 ). This is the discrepancy between the measurements 
noted above; a possible resolution is discussed in Sec. IIV Al 

As a test of the sensitivity of X^sr to the data fitting procedure, the functional form 
G(t) cos(ut + 0), with a damped Gaussian envelope 

G(t) = e- m exp(-|(T 2 t 2 ), (10) 

was also used to fit the TF-/1SR data from both the normal and superconducting states. 
In the normal state, where this function provides a good fit, the Gaussian relaxation arises 
from quasistatic Sb nuclear dipolar fields. The exponential damping increases with increasing 
field and is mainly due to normal-state susceptibility inhomogeneity; dynamic fluctuations 
of hyperfine-enhanced 141 Pr nuclear spins 48 play a small role at low fields. Below T c the 
exponential rate W was fixed at the normal-state value for each field, so that the temperature 
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dependence of the Gaussian rate a reflects the effect of the superconducting state. Some 
such procedure is necessitated by the strong statistical correlation between W and a in 
Eq. (1101) ; the time constant and the shape of the relaxation function are influenced by both 
of these parameters, so that correlations between them can result from small systematic 
errors. The principal justification for this ad hoc fixing of W is the insensitivity of the 
superconducting-state results to details of the fitting function [Eq. (flQj) ] discussed below. 

Determination of the vortex-state field distribution width requires correction for the 
normal-state relaxation. We take the superconducting-state Gaussian rate a s to be given by 
o 2 = a 2 — a 2 , where a n is the normal-state rate.™ 

Equation ([3]) relates the second moment SB 2 of field distribution to A in the London limit. 
The second moment of the corresponding fi + frequency distribution is 5uj 2 = / y 2 5B 2 , where 
7 M is the /i + gyromagnetic ratio. Then the estimated penetration depth A from Eq. ([3]) is 



A (/im) = 0.328/v^u; (/xs" 1 ). (11) 

As noted above the rms width o~ s of the best-fit Gaussian is not necessarily 5uj, so that 
replacement of 5u in Eq. (ITT]) by o~ s is not completely justified. Nevertheless a s should 
scale with 5uj, and within its range of validity Eq. fTTTj) should give the correct temperature 
dependence of A. This is because under these circumstances effects of nonzero £ are restricted 
to the high-field tail of field distribution, which is not heavily weighted in a Gaussian fit 



(cf. Fig. 1 of Ref. |50|). PrOs 4 Sbi 2 is a strongly type-II superconductor (GL k = A/£ ~ 30, 



Refs. 



10 and 



251 ) . and this picture should be applicable. The temperature dependencies of 
A^sr obtained from Gaussian fits and Eq. fTTTj) (circles) and from GL model fits (squares)^, 
are shown in the inset of Fig. [T] for an applied field of 10 mT. The results agree very well, 
giving strong evidence that X^sr is robust with respect to very different fitting strategies. 

Figure [2] gives the temperature dependence of the corrected superconducting-state fi + 
spin relaxation rate o~ s (T), normalized to the zero-temperature value cr(0), for fioH = 10, 
20, and 100 mT. The relaxation rates are well described by a power-law temperature 
dependence cr(T) = er(0)[l — (T/T c ) n ] (curves in Fig. [2]) over the entire temperature range 
T < T C (H), with the parameter values given in Table [B For all of the present data n > 4, 
and for fioH = 10 and 20 mT o~ s {T) is nearly temperature-independent below ~0.8 K. These 
features suggest the absence of gap nodes^^ 1 ^ at these fields. We should note, however, that 
such power-law fits have no clear physical interpretation and, furthermore, are dominated 
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4 12 




T/T C {H) 



FIG. 2: (color online) Normalized corrected superconducting-state /i + spin relaxation 
rates a s (T)/a s (0) in PrOs4St>i2. (a): applied field fj>oH = 10 mT (circles), (b): hqH = 20 mT 
(squares), (c): {jlqH = 100 mT (open triangles). Solid triangles: 100 mT data from Ref. !38. Curves: 
fits of the power law a(T)/a(0) = [1 — [T/T c ) n ] to the TF-/xSR data. Diamonds: inverse-square 



Meissner-state penetration depth A 2 (0)/A 2 (T) from rf data (Ref. 



m. 



TABLE I: Parameters from power-law fits cr s (T) = o" s (0)[l — (T/T c ) n ] to the corrected 
superconducting-state fi + spin relaxation rate cr s (T) in PrOs4St>i2. 



Field (mT) 



<T 8 (0) (MS" 1 ) 



Tc(H) (K) 



Exponent n 



10 

20 
100 
100" 



0.82(1) 
0.77(1) 
0.80(2) 
0.75(2) 



1.83(1) 
1.79(2) 
1.68(4) 
1.80(2) 



4.5(3) 
4.3(4) 
4.1(6) 
2.9(3) 



"From Ref. 



38. 



by the behavior of cr s (T) at higher temperatures, which is not the concern of this article. 
The large values of n are merely signatures of the near temperature independence at low 
temperatures. 

The normalized inverse-square penetration depth A 2 (0)/A 2 (T) from rf measurements us- 
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ing A(0) = 0.353 jttm (Fig. [TJ), which is equal to a s (T) / a s (Q) in London superconductor as 
noted above, is also shown in Fig. [2] for comparison. For low fields a small but clear dis- 
crepancy at low temperatures can be seen: the relaxation rates decrease significantly less 
rapidly than A 2 (0)/A 2 (T) with increasing temperature. This is of course the same discrep- 
ancy shown in Fig. [TJ As noted above the effective penetration depth from TF-/zSR can 
be modified by vortex interactions at higher fields;— the decrease of n with increasing field 
(Tabled]) may reflect such an effect. 



38 taken 



Figure W[c) compares our TF-/1SR data for /j, H = 100 mT with those from Ref. 
at the same field. There is reasonable agreement between all the data at low temperatures, 
but the exponential damping [Eq. (ITU]) ] increases with field and begins to dominate the 
relaxation, thereby increasing the error in the Gaussian relaxation rate. The resultant 
scatter in the /iSR data leads to noticeable differences in fit-parameter values between the 
two experiments (Table [I]), so that the situation for fi H = 100 mT is uncertain. 

These results, together with a previous TF-/zSR study on a different sample,— indicate a 
generally flatter temperature dependence of \^,sr{T) at low temperatures compared to that 
of A r f(T). As a whole the data suggest that the A M sR-A r f discrepancy is intrinsic and robust 
for low fields (< 20 mT), where the /zSR penetration depth determination is most reliable.— 

It can also be seen in Table [I] that the zero-temperature relaxation rate is essentially 
independent of field. In an isotropic superconductor such as cubic PrOs4Sbi2 vortex-lattice 
disorder is expected to increase the low-field rate; increasing field (increasing vortex density) 
then decreases the rate as intervortex interactions stabilize the lattice.— ^2 Thus the field 
independence of the low-temperature rate indicates a substantially ordered vortex lattice at 
all fields, in which case the temperature dependence of the rate is controlled solely by the 
temperature dependence of the effective penetration depth. 



B. Magnetic susceptibility 

The temperature dependence of the normal-state magnetic susceptibility of the PrOs 4 Sb 12 
sample used in the TF-/zSR experiments has been determined using a commercial SQUID 
magnetometer. Figure [3] shows the measured bulk susceptibility (circles). The data were fit 
to the model of CEF-split Pr 3+ ions with cubic Oh symmetry interacting via a molecular 
field as discussed in Sec. II CI The resulting temperature dependence is given by the curve in 
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1 — I 




Temperature (K) 

FIG. 3: Temperature dependence of measured bulk susceptibility (circles) and best fit to model of 
interacting CEF-split Pr 3+ ions (curve) in PrOs4Sbi2- 

Fig. [31 The fit values of the CEF parameters in Eq. ([7j) are listed in Table [TTJ The fit value 



TABLE II: Fit values of CEF parameters in meV for PrOs4Sbi2- Oh- cubic symmetry (no \jl 
perturbation). Ci v : tetragonal symmetry (^ + perturbation, cf. Sec. II C|) . 



o h 


R° 

Civ 


R° 

o h 


Civ 


Bi 

Oh Ci V 


R° 

o h 


Civ 


BI 

o h 


Civ 





-0.18 


0.0154 


0.0153 


0.0771 0.0771 


0.0007 


0.0007 


-0.0154 


-0.0154 



of the molecular field parameter i was found to be —1.67 mole/emu, which is close to the 
value of —2.54 mole/emu found by Tayama et al. 51 

C. Normal-state muon Knight shift 

We have performed TF-/1SR experiments at applied field n$H = 1.0 T in the normal 
state of PrOs4Sbi2 . Since the fi + frequency shift is proportional to the magnetic field at the 
fi + site, it is resolved better at higher fields. Fig. @|a) shows the temperature dependence 
of the muon-spin precession frequencies from the sample and an Ag reference. Silver was 
used because Ag nuclear dipole fields are weak, leading to a well-defined reference frequency 
(narrow line), and the /i + Knight shift in Ag is known (94 ppm).^ One can extract the /i + 
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frequency to an accuracy of ~100 ppm.— The temperature dependence of the relative 
frequency shift K* is shown in Fig. H] (b). 



PrOs „Sb u H=1.0T 

4 12 H) 




100 200 

Temperature T(K) 



FIG. 4: (a): Temperature dependencies of // + precession frequencies in PrOs4St>i2 (uj and an 
silver (Ag) reference (u Te {). (b): Temperature dependence of /i + frequency shift K* in the normal 
state, (c): Temperature dependence of relative rms linewidth 5K^ = cr/uj Te f, where a is the /i + 
Gaussian relaxation rate. 

Figure [5] (triangles) gives the dependence of K* on the bulk molar susceptibility Xmoi k m 
PrOs4Sbi2 in the normal state, with temperature an implicit variable. K* and x^t were 
measured in the same sample. As discussed in Sect. II Bl K* should be corrected for the effect 
of the Lorentz and demagnetization fields 4.7r(l/3 — D)M, where M is the magnetization. 
The demagnetization factor D was estimated from (i) the ratio of the height and width 
of the samples,— with the magnetic field applied perpendicular to the flat faces of the 
rectangular crystallites, and (ii) the ratio / of individual crystal volume to the sample as a 
whole, taking spaces between the crystallites into account. We estimate the demagnetization 
factor for the entire sample D samp = 0.824, the demagnetization factor for individual crystals 
D CTys = 0.365,^ and / = 0.95. Then the effective value is D — fD 8amp +(l-D)n ciys = 0.80.^ 
Therefore, K DL = 4tt(1/3 - n)xy u lk = -4vr(0.47)xy lk , where Xy lk is the bulk susceptibility 
per unit volume, or K DL = — 0.0243Xmoi k - We subtract this from K* to obtain the corrected 
/x + Knight shift (Fig. O circles). A linear relation is obtained for x^, lk ^ 8-0 x 10 -2 

15 




FIG. 5: (color online) Dependence of \x + frequency shift K* (triangles) and corrected [i + Knight 
shift (circles) on bulk molar susceptibility X^T m PrOs4Sbi2, applied field ^qH = 1 T. The 
straight line K = -Kbulk ( see text) is a fit to the data in the region 6.3 K < T < 250 K. 

emu mole -1 (T > 6.3 K). This implies that at high temperatures the fi + shift samples the 
same electrons that produce the large temperature-dependent bulk susceptibility component. 
Such behavior is generally expected from the "bulk" Knight shift i^buik m the absence of 
local perturbations. Only a small deviation from linearity (a single point) appears below 
6.3 K. 

As discussed in Sec. II B\ the average shift arises solely from the hyperfme contact field at 
the interstitial /i + for a cubic lattice. The dipole-dipole interaction, however, may split or 
broaden the /z + line. Assuming the most probable /i + stopping site (|, 0, 0.15) determined by 
Aoki et a/.,—, we obtain a calculated dipolar Knight shift tensor with the following principal 
values for field directions (x,y,z) (cf. first column of Table HTT|) : K%L = —2.70 x 10 -2 Xmoh 
K y dip = 1.121 x 10- 2 Xmoi, and K z dip = 1.58 x 10 -2 Xmoi- Thus in the absence of additional 
broadening the fi + spectrum should be split into three lines of equal weight. The observed 
relative rms linewidth 5K^ = 4.68 x 10~ 2 Xmoi (data of Fig. BJc); dependence on Xmoi not 
shown) is, however, larger than the rms spread Ai-Qip = \j[i.K dip ) 2 + {K dip ) 2 + {K dip ) 2 }/3 = 
1.92 x 10~ 2 Xmoi- Thus the dipolar splitting cannot be resolved, although it contributes 
significantly to the linewidth. 
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IV. DISCUSSION 



A. Multiband Superconductivity 

We first consider the discrepancy in measured magnetic penetration depth between the 
TF-/iSR and rf experiments. Recently, extreme MBSC was found in PrOs4Sbi2 from heat 
transport measurements by Seyfarth et al.- 1 ^ Their thermal conductivity and other data are 
explained by small and large gaps A s , A; on different sheets of the Fermi surface, together 
with different Fermi velocities vp s , vpi and coherence lengths £ s j ~ Hvp S) i/A s j. A crossover 
field H^2, which corresponds to the overlap of the vortex core electronic structure due to the 
small-gap band, is given by H s c2 = $ /27r£| ~ 10 mT,- which is of the order of the lower 
critical field H c \. Microscopic calculations of the local DOS in a two-band superconductor^ 
suggest that the small-gap band induces spatially extended QP states at low field. For 
H > H^2, these loosely bound states overlap and become delocalized, with the dominant 
contribution to the DOS coming from the large gap band. In PrOs4Sbi2 this high-field region 
includes most of the vortex state; it is in this sense that PrOs4Sbi2 is an extreme multiband 
superconductor. 

Our TF-/iSR measurements were performed for applied field [IqH > 10 mT « H^ 2 - Then 
the small-gap states and their contributions to screening supercurrents are nearly uniform, 
the vortex-state field inhomogeneity is mainly due to large-gap suppercurrents, and A ex- 
hibits an activated temperature dependence if the large gap is nodeless. In contrast, the rf 
measurement of the surface penetration depth was performed in the Meissner state. Both 
large- and small-gap Cooper pairs contribute to the superfluid density. Its temperature 
dependence is controlled by both small- and large-gap superfluid densities; the small-gap 
contribution dominates the temperature dependence at low temperatures. Then the tem- 
perature dependence of the penetration depth in the vortex state from TF-/iSR is weaker 
than that in the Meissner state from rf measurements. It should be noted that in this sce- 
nario the TF-/iSR measurements give no information on the nodal structure of the small 
gap. 
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B. Low-temperature Knight shift 

In Fig. [5j a small deviation from the linear relation of the /x + Knight shift vs bulk 
paramagnetic susceptibility ^ bulk appears below 6.3 K in PrOs4Sb 12 . Kaplan, Schenck et 
al. have reported deviations from linear K-\ relations from TF-/1SR measurements of the 
/i + Knight shift in single crystals of PrNig^^ 1 ^ and Prln3.— These were attributed to 
induced changes of the low-temperature susceptibility due to modification of the CEF of 
neighboring Pr 3+ ions. We argue that the /z + 's perturbing effect is small in PrOs4Sbi2 for 
the following reasons: 

• The deviation of K^{x) from linearity in PrOs 4 Sb 12 is very small (Fig. E]). 

• The superconducting transition temperature measured by TF-/iSR is consistent with 
the bulk superconducting value (see Fig. |2J), so that this signature of superconductivity 
is not affected by the fi + charge. 

• Pr 3+ ions are considerably more dilute in PrOs4Sbi2 (2.48 x 10 21 cm -3 ) than in PrNis 
(11. 8x 10 21 cm -3 ) or Prln3 (9.81 x 10 21 cm -3 ); the nearest-neighbor /i + -Pr 3+ distance is 
about 2 A longer in PrOs4Sbi2 than in PrNi 5 or Prln 3 . Thus for comparable screening 
lengths the /i + electric field at neighboring Pr 3+ sites in PrOs4Sbi2 is more completely 
screened by the conduction electrons. Similarly, the experimental slope dKbuik/ dxmoi = 
—4.8 x 10 -2 mol cm -3 in PrOs4Sbi2 (Fig. [5]) is smaller than the corresponding values in 
PrNis and Prln3 by more than an order of magnitude, indicating the RKKY coupling 
between fi + spins and neighboring Pr 3+ ions is correspondingly weaker in PrOs4Sbi2- 

• Schenck et al— concluded from ZF-/iSR studies that the spin fluctuations of the Pr 3+ 
electronic moments in PrNi 5 are slowed down and exhibit quasistatic behavior at low 
temperatures. They suggested that this behavior is related to /^-induced modification 
of the low temperature susceptibility of the neighboring Pr ions. We did not observe 
such large widths of the magnetic field distribution at low temperatures [see Fig.|U(c)]. 
Our ZF- and LF-/1SR measurements^ suggest instead that the observed dynamic /i + 
relaxation is due to the 141 Pr nuclear spin system with an hyperfine-enhanced effective 
nuclear moment. 

We now discuss the small deviation from the linear K-\ relation that appears in Fig. [5] 
below 6.3 K. It should be noted that muon-induced CEF modification is not the only mech- 
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anism for a nonlinear K(x). It has been recognized for some time that such behavior, first 
observed in NMR Knight shifts, 56 might arise from effects such as (1) temperature-dependent 
modification of the hyperfine coupling or (2) the onset of a new susceptibility component at 
low temperatures. Thermal depopulation of excited CEF states at low temperatures, with 
consequent modification of the coupling,— is an example of a Type-(l) mechanism. Non- 
linear K(x) behavior in heavy-fermion materials has recently been ascribed 5 ^ to a Type- (2) 
effect, related to a two-fluid description of the heavy-fermion state.— Our concern here is 
the possible influence of the muon charge on our vortex-state TF-/zSR results in PrOs 4 Sbi 2 ; 
if other mechanisms were responsible for the small K-\ nonlinearity there would be no 
evidence for such influence. 

We therefore suppose the nonlinearity in PrOs4Sbi2 is due solely to the muon's perturbing 
effect, i.e., the local modification of the Pr 3+ CEF splitting changes the local Pr 3+ suscep- 
tibility, so that the fi + Knight shift is no longer monitoring the bulk susceptibility. Thus 
we can use the /i + Knight shift data to determine this change. To obtain a quantitative 
determination, we use the approach of Schenck et a/.,— >2I making the following assumptions: 
(i) only the susceptibilities of the two nearest Pr 3+ ions are affected by the muon, and (ii) 
only the two nearest Pr 3+ ions contribute to the contact interaction at the fi + site. Based 
on these assumptions, with the principal axes of x chosen as the coordinate frame (x, y, z), 
the Knight shift with the external field in the i direction (i — x,y or z) may be written as 

K % = (A" ip)NN + Ax> n )xj ocal + ^di P ,(i-NN)X bulk - (12) 

Here xl° cal is the altered susceptibility of the two nearest Pr 3+ ions, and ^ bulk is the unper- 
turbed bulk Pr susceptibility. The subscript NN signifies that the sum in only includes 
the two nearest Pr neighbors, and the subscript 1 — NN indicates summation over all Pr 3+ 
ions in the Lorentz sphere other than the nearest neighbors. For the most probable fi + 



site (~, 0, 0.15) (Ref. Il7|), the calculated values of A^i P , Aiip.NN, and A^i-nn are listed in 
Table IIHI We take A con to be the slope of K^-\ in Fig. [51 since the tensor A dip is traceless 
and the sum £\ A% ip vanishes as can be seen from Table IIHI 

In order to obtain the values of Xi° cal , we must modify Hcef to produce the required level 
changes. In the presence of a muon, the symmetry of the neighboring Pr 3+ ions is not the 
original cubic symmetry, but rather a very approximate tetragonal (C& v ) symmetry. The 
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TABLE III: Calculated values of /x + -Pr 3+ coupling constants A^p, Adi P: i-NN, and A COIl (see text) 
in PrOs4St>i2. 

i A% ip (mol/cm 3 ) ^dip.NN (mol/cm 3 ) A % p ,i-nn ( m ol/cm 3 ) A con (mol/cm 3 ) 

-5.08 x 2.37xl0- 2 

2.90xl0~ 2 -1.78xl0~ 2 4.8xl0" 2 

2.18xl0~ 2 -0.60 x 10" 2 



general Hamiltonian describing this tetragonal symmetry is 34 ' 35 

^ C EF,tet = BlOl + BlOl + B\0\ + B%0% + BlOl (13) 

which has five free parameters, compared to two in the original cubic Hamiltonian. In 
principle, the muon may modify the charge distribution around the Pr ions in every direction. 
However, if we only consider changes occurring along the Pr 3+ -/i + axis for simplicity, it has 
been shown that B\ and B\ remain at their cubic values.— This reduces the number of the 
free parameters from five to three. Hence, using Eq. (fl~2l . we can fit the /i + Knight shift 
by choosing appropriate values of £> 2 , B\, and Bq. Our theoretical fits in PrOs4Sbi 2 are 
shown in Fig. El The anisotropy is consistent with the new CEF and Eq. ( fl2]) . decreases with 
increasing temperature, and is well within the spread SK^ of /i + shifts (shown as "error bars" 
in Fig. [6]). The data are very well described by the average shift K w = (K x + K y + K z )/3, 
which is almost the same as i^buik from Fig. [5] that assumes no fi + perturbing effect. The 
fitted CEF parameters are listed in Table [TT] under the columns with heading C<± v . The 
changes in the parameters, which are very small and probably not statistically significant, 
give rise to a correspondingly small rearrangement of the local energy levels. A comparison 
of bulk and locally perturbed energy level scheme of PrOs4Sb 12 is shown in Fig. [71 We see 
that the change of splitting between the ground state and first excited state energies is only 
0.3 K, and conclude that the perturbing effect of the /i + charge in PrOs 4 Sbi 2 is negligibly 
small. 



x -2.70 x 10~ 2 

y 1.121xl0~ 2 
z 1.58xl0" 2 



V. CONCLUSIONS 



In this paper we have presented TF-/iSR results in the first Pr-based heavy- fermion 
superconductor PrOs 4 Sbi 2 . It is found experimentally that the effective penetration depth 
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FIG. 6: (color online) Temperature dependence of /x + Knight shift in PrOs4Sbi2- Circles: fi + 
Knight shift K„ from fig. "Error bars:" Gaussian linewidth 5K^ [Fig. H(c)]. Dash-dot curve: 
calculated K x . Dashed curve: calculated K y . Dot curve: calculated K z . Solid curve: average 
shift -K" av = (K x + K y + K z )/3. Dash-dot-dot curve: Kbulk from linear fit (Fig. [5]), expected for no 
/i + perturbing effect. 
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FIG. 7: Original bulk energy-level scheme of the Pr ions in PrOs4Sbi2, compared to the local level 
scheme induced by 

A e ff is temperature-independent in the vortex state for low temperatures, consistent with 
a nonzero gap for quasiparticle excitations. In contrast, the temperature dependence of 
the penetration depth from rf inductive measurements suggests point nodes in the gap. 
This discrepancy can be resolved in a scenario based on the recent discovery of two-band 



5.5 K 
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superconductivity in PrOs4Sbi 2 . 

The temperature dependence of the normal-state fj, + Knight shift in PrOs 4 Sbi 2 reveals 
a linear scaling of the Knight shift with the bulk magnetic susceptibility at high tempera- 
tures. A very small deviation from the linear relation appears in PrOs4Sb 12 below 6.3 K. 
Such deviations have been explained by fi + induced modifications of the susceptibility of 
neighboring Pr 3+ ions due to a change of the CEF splitting. Our data indicate, however, 
that this modification is very small in PrOs4Sbi 2 . A model calculation based on CEF the- 
ory and the associated perturbed electronic energy levels confirms the smallness of the /z + 
perturbation effect. Therefore it is unlikely that the discrepancy between TF-/iSR and rf 
inductive measurements is caused by the /i + charge or, more generally, that the /z + charge 
influences the role of the Pr 3+ CEF states in the superconductivity of PrOs 4 Sb 12 . 
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